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We present the results of a weakly modeled burst search for gravitational waves from mergers of 
non-spinning intermediate mass black holes (IMBH) in the total mass range 10CM50 M© and with 
the component mass ratios between 1:1 and 4:1. The search was conducted on data collected by the 
LIGO and Virgo detectors between November of 2005 and October of 2007. No plausible signals 
were observed by the search which constrains the astrophysical rates of the IMBH mergers as a 
function of the component masses. In the most efficiently detected bin centered on 88 + 88 Mq , for 
non-spinning sources, the rate density upper limit is 0.13 per Mpc 3 per Myr at the 90% confidence 
level. 

PACS numbers: 95.85.Sz, 04.80.Nn 


I. INTRODUCTION 

Emission of gravitational waves (GW) via strong gen- 
eral relativistic processes between two compact objects 
(black holes and/or neutron stars) is the hallmark of com- 
pact binary coalescence (CBC). Binary black holes, a par- 


ticular class of CBC sources, have been one of the main 
detection targets of ground based gravitational wave de- 
tectors since the inception of large wide-band interferom- 
eters [1-3], This paper presents the results of a search 
for gravitational waves from the coalescence of interme- 
diate mass black holes (IMBH). The search used data 
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collected by the Laser Interferometer Gravitational- Wave 
Observatory (LIGO) during its fifth science run (S5) from 
November 2005 to October 2007 [1] and by the Virgo GW 
interferometer [4] , which commenced its first science run 
(VSR1) in May 2007 and operated in coincidence with 
LIGO. 

The coalescence of compact binaries is generally di- 
vided into three stages: the inspiral, merger, and ring- 
down. Gravitational waves from the inspiral stage are 
quasi-periodic “chirp” signals of increasing frequency and 
amplitude which are well described by analytical post- 
Newtonian (PN) models [5-8] before the binary evolution 
reaches the inner-most stable circular orbit (ISCO). Near 
the ISCO, the strong gravitational interaction no longer 
allows for a stable orbit and the two black holes merge to- 
gether to form a single black hole. After the merger stage, 
the newly born perturbed black hole emits gravitational 
waves via exponentially damped quasi-normal modes in 
the ringdown stage. The merger and ringdown stages 
of the GW signal are important for detection of IMBH 
sources because the characteristic frequencies of the in- 
spiral stage are usually outside of the sensitivity band 
of ground-based GW interferometers. Recent progress in 
numerical relativity (NR) has expanded the understand- 
ing of binary black hole systems through the merger and 
ringdown stages [9-14] allowing calculation of the full 
inspiral-merger-ringdown (IMR) waveforms. 

Several matched-filter searches have been conducted 
for CBC sources consisting of total masses less than 
35 M© [15-17] with inspiral templates. One more search 
in the total mass range 25-100 M© (where the contribu- 
tion of the inspiral stage is still dominant) was performed 
with IMR templates [18]. In order to identify GW events 
in the noisy data, these searches rely on the generation 
of template banks from the signal model. Currently, the 
generation of complete and accurate template banks in 
the mass region above 100 M© is challenging. Therefore, 
for the IMBH search reported here we used the Coher- 
ent WaveBurst algorithm [19-21] which is designed for 
detection of un-modeled burst signals and does not re- 
quire a priori knowledge of the signal waveforms. How- 
ever, due to the lack of model constraints, generic burst 
searches are usually more affected by the background 
than matched-filter searches. To improve the rejection 
of background events, the CWB algorithm can enforce a 
constraint on the waveform polarization [22] . Such a con- 
strained burst algorithm can be used to search for IMBH 
coalescences without the need of template banks while 
still achieving nearly the same detection sensitivity. 


A. Intermediate Mass Black Hole Formation 

IMBHs have been posited to complete the black hole 
mass hierarchy. As such, IMBHs cover several decades in 
the black hole mass spectrum between stellar mass black 
holes of a few tens of M © , formed from star collapse, and 
super-massive black holes of 10 5 M© or more present in 


the center of galaxies. Some models of IMBH formation 
include runaway stellar collision scenarios [23] in globu- 
lar clusters (GC). One model proposes that lower mass 
single IMBHs could be formed by the stalled supernova 
of early Population III stars [24, 25]. Another model [26] 
studies the progressive accumulation of mass into a large 
(> 50 M© ) seed black hole via coalescence of a popu- 
lation of smaller black holes. However, the existence of 
binaries with IMBH components remains uncertain since 
stellar winds may stall the growth of the IMBH progeni- 
tors in the runaway collision scenario [27], or the merger 
recoil may also eject a newly formed black hole from the 
globular cluster [28, 29]. 

IMBHs have been searched for via conventional as- 
tronomy, and a few candidates exist [30]. It has been 
suggested that IMBHs are the engines powering ultra- 
luminous X-ray (ULX) sources [31, 32] such as M82 X- 
1 [33, 34] or NGC 1313 X-2 [35]. Most models agree 
that the primary hosts of these objects would be globu- 
lar clusters [36-38]. These objects are thought to grow 
from accretion of smaller compact objects [39] , and there- 
fore, IMBHs could be a prime candidate for the detection 
of GW by the coalescence of solar mass objects into the 
central BH. The detection of an IMBH binary would not 
just represent the first detection of GW, but could have 
important consequences for theories about the formation 
of super-massive black holes and the dynamics and evo- 
lution of globular clusters [38, 40]. 

In this search, we focus on the IMBH binary systems 
with total masses between 100-450 M©and component 
black hole masses with the mass ratios between 1:1 and 
4:1. The expected GW emission from these sources is in 
the frequency band between tens and few hundred Hz. 
This frequency band includes the most sensitive band of 
the initial ground based GW detectors. Those IMBH 
systems considered in this search contain most of the de- 
tectable signal power in the merger and ringdown because 
the power emitted during the earlier inspiral stage is in 
the frequency band below 40 Hz where there is a rapid 
deterioration of the LIGO/Virgo network sensitivity due 
to seismic noise. Above 450 M© , the power emitted is 
no longer present at accessible frequencies. 

Concerning the rate of IMBH-IMBH coales- 
cence, the upper limit has been estimated at 0.07 
GC _1 Gyr _1 [41]. Using an astrophysical source 
density of 0.3 GC Mpc -3 [42], this corresponds to 
2xl0“ 5 Mpc -3 Myr _1 . If intermediate mass ratio 
(mass ratios of 10:1 or greater) inspirals onto IMBH 
are also considered, the rate is conceivably as high as 
3 GC _1 Gyr -1 (9xl0“ 4 Mpc^Myr” 1 ). The detection 
rate estimates for the IMBH systems considered in 
this search are much smaller than 1 yr” 1 . However, 
predicted detection rates for second generation detectors 
such as Advanced LIGO and Virgo increase by orders 
of magnitude over initial detectors as their proposed 
designs include better sensitivity at comparatively lower 
frequencies. 
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II. EXPERIMENT 

Five GW detectors were operating during the 
S5/VSR1 runs: two detectors (4 km detector HI and 
2 km detector H2) at the LIGO site in Hanford, Wash- 
ington, another 4 km LIGO detector (LI) in Livingston, 
Louisiana, the 3 km Virgo detector (VI) in Cascina, Italy, 
and the 600 m GEO600 detector in Hannover, Germany. 
The GEO600 detector had a significantly lower sensitiv- 
ity to the IMBH sources than the other four detectors 
and therefore it was not considered in this search. Due to 
limited detector duty cycles, there were several network 
configurations consisting of two to four detectors oper- 
ating in coincidence. In this search we considered two 
networks with the most accumulated observation time: 
the three-fold network L1IIIII2 and the four-fold network 
L1H1H2V1. 

Not all data which was collected by the detectors is 
used in the analysis. Extensive studies [43, 44] have been 
performed to identify (flag) data segments with high seis- 
mic activity, large mechanical disturbances, and a high 
rate of environmental and instrumental transients. These 
data quality flags are nearly identical to those used for 
the S5/VSR1 all-sky burst analysis [20, 21]. Data quality 
flags are classified into different categories, starting with 
the initial flags selecting data segments used by the search 
algorithm. Further data quality flags are imposed on all 
events emerging from the search algorithm, including a 
set of event vetoes derived from well known correlations 
between the GW data channel and the auxiliary chan- 
nels. All events passing these checks are considered as 
detection candidates. Finally, a set of data quality flags 
is used to remove events with weaker environmental and 
instrumental correlations. The set of events passing the 
final checks is then used for estimation of the astrophys- 
ical rate limits. Table I shows the total observation time 
for the network configurations used in the search after all 
data quality flags are applied. 

TABLE I: Summary of each network’s analyzed observation 
time after all data quality flags are applied. 

detector network observation time (yr) 

L1H1H2V1 0.16 

L1H1H2 0.65 


III. INTERMEDIATE MASS BINARY BLACK 
HOLE SEARCH 

A. Search Algorithm 

The IMBH search is based on the Coherent Wave- 
Burst (CWB) algorithm [45] which has been used in the 
S5/VSR1 burst searches [20, 21]. The CWB algorithm 
performs a constrained likelihood analysis [46] of the 


network data stream, reconstructing detector responses 
to an anticipated GW signal. The residual data (null 
stream), obtained after subtraction of the estimated de- 
tector responses from the original data, represents the 
reconstructed network noise. Along with the reconstruc- 
tion of un-modeled burst signals, which imply random 
polarization, CWB performs likelihood analysis of sig- 
nals with other polarization states, including elliptical, 
linear, and circular polarizations. 

In this search, we use the elliptical polarization con- 
straint. The details of the likelihood analysis and the 
elliptical polarization constraint are presented in Ap- 
pendix A. Though not completely generic, the constraint 
improves rejection of those background events originat- 
ing from the random coincidence of the environmental 
and instrumental transients in the detectors. In general, 
the signal polarization may evolve as a function of time. 
For example, spinning black hole systems have slowly 
evolving, large polarization changes that track the pre- 
cession of the orbital angular momentum [47] . Even non- 
spinning black hole systems have small, rapidly oscillat- 
ing polarization changes as different multipolar orders in- 
terfere constructively and destructively [48] . Either case 
will introduce a time dependence on the waveform polar- 
ization. However, the effects of the spin-orbit coupling 
become significant only for a subset of black hole sys- 
tems with spinning components. Even for these signals, 
a significant fraction of the band-limited power can be 
associated with some instantaneous polarization. Thus, 
the constraint should not significantly affect the detec- 
tion efficiency of IMBH sources. 

Three major statistics — obtained as the result of 
the likelihood analysis — are used for selection of re- 
constructed events: the network correlation coefficient 
cc, the network energy disbalance A ne t and the coher- 
ent network amplitude rj (see Appendix A). The statis- 
tics cc and A net are used to characterize the conformance 
of identified events with the signal model and its con- 
straints. A low value of cc < 1 is typical for background 
events which tend to have a large residual energy and a 
small coherent energy of reconstructed signals. On the 
contrary, a genuine GW event is characterized by a value 
of cc close to unity. The energy disbalance A ne t identi- 
fies the un-physical solutions of the likelihood functional 
which are typical for spurious events. A significant devi- 
ation of A ne t from zero is an indication that the energy 
of the reconstructed response is significantly larger than 
the energy of the data stream in at least one detector. 
Table II shows the cc and A net thresholds used in the 
analysis. 


TABLE II: Post-production selection cuts: candidate events 
are selected if A nc t and cc are, respectively, less and greater 
than indicated thresholds. 


network 

H1H2L1V1 H1H2L1 

dual stream energy disbalance ( A ne t) 
network correlation coefficient (cc) 

0.2 

0.6 

0.15 

0.70 
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The coherent network amplitude r? is the main CWB 
detection statistic. It is proportional to the signal-to- 
noise ratio (SNR) and is used to rank selected events 
and e&tablish their significance against a sample of back- 
ground events. 

B. Background Estimation 

We estimate the false alarm rate of events originat- 
ing from the detector noise by introducing artificial time 
shifts (far exceeding the intersite light travel time) be- 
tween the data from different sites before using the search 
algorithm. This procedure assumes that the noise in- 
duced events are not correlated between the sites. Events 
obtained from the time-shifted data represent the search 
background sample. Data from different detectors is 
shifted by integer multiples of one second per time-shift 
configuration (time lag). The H1H2L1 network had a 
total of 600 time lags performed, and the H1H2L1V1 
network had a total of 1000 time lags including the fore- 
ground (zero lag) configuration. In total, this procedure 
accumulated 569 years of effective background live time 
for the three detector network and 180 years of back- 
ground live time for the four detector network. The 
background events that survived the data quality and 
the analysis selection cuts (see Table II) are used for cal- 
culation of the significance of candidate events and the 
false alarm density statistic described in section IV C. 


IV. SIMULATIONS 

To characterize the detection efficiency of the search in 
the parameter space of potential IMBH sources, extensive 
simulation studies were performed with different families 
of the IMR waveforms. These studies were made to de- 
termine a sensitivity volume of the search, also called 
visible volume, assuming that the IMBH sources are dis- 
tributed uniformly in space. In order to calculate the 
visible volume a Monte Carlo detection efficiency study 
was performed by adding into the data waveforms drawn 
randomly from the physical parameter space which we 
consider. The simulated detector responses were injected 
via software into detector data and the search algorithm 
was used to identify the injections. A large sample of 
waveforms for each network configuration was generated 
to sufficiently cover the parameter space of the IMBH 
sources presented in Table III. The simulated waveforms 
were distributed in a spherical volume with a radius of 2 
Gpc and a uniform distribution over the source inclina- 
tion and polarization angles, and sky locations. 

In this simulation, redshift corrections were neglected 
because very few injections are placed (and detected) at 
the distances which would require consideration of this 
effect. Spin of the component black holes was not con- 
sidered as well, but a discussion of potential spin effects 
is presented in section VI. 


TABLE III: Summary of injected waveform parameters. 
Total Mass (M© ) 100 - 450 
Mass Ratio 1-4 

Distance (Mpc) 0 - 2000 


A. Simulated Waveforms 

Most of the previous template based searches [16, 17, 
20, 49] used only inspiral or ringdown templates [50] to 
do simulation studies. As the total mass of the system 
increases, the analytical PN inspiral waveforms become 
inadequate because only the merger and ringdown waves 
have significant power in the sensitive band of the detec- 
tors. For this reason, this search (and a template search 
for binary black holes with the total mass between 25 
and 100 M 0 [18]) uses the full IMR waveforms from 
two different families: the Effective One Body Numer- 
ical Relativity (EOBNR) family [9, 10, 12, 51] and the 
IMRPhenom family [11, 52, 53]. The EOBNR waveform 
family uses the Effective One Body (EOB) Hamiltonian 
to evolve the binary system up to the merger. The EOB 
approach is able to simulate the dynamics of the plunge 
into merger of the binary black hole system through 3 PN 
order. Further accuracy has been obtained by the use of 
“pseudo” 4 PN terms motivated from the results of nu- 
merical relativity simulations. To complete the evolution 
from the plunge-merger to ringdown, a superposition of 
the ringdown frequency modes is matched to the end of 
the merger. The IMRPhenom family is constructed by 
matching 3.5 PN order analytical inspirai waveforms to 
the corresponding NR merger waveforms to make “hy- 
brid” waveforms. These hybrid waveforms are then ex- 
trapolated to form a full waveform family in the Fourier 
domain. In contrast to [18] which constructed template 
banks from the EOBNR family but used both families 
for detection efficiency studies, this search uses only the 
EOBNR family for efficiency studies. While not used for 
a detailed simulation, the IMRPhenom family was used 
in the IMBH analysis to cross-check the validity of results 
obtained with the EOBNR family. 

B. Visible Volume 

In general, the visible volume [18] is a function of the 
component masses (mi, m 2 ) of the binary system. It can 
be calculated as 

/»oo 

Vvi S (mi,m 2 ,J?) = 4 tt / e(r,mi,m 2 ,??)r 2 dr . (4.1) 

Jo 

where e is the detection efficiency of the search, which 
is also a function of the distance to the source r. The 
visible volume is calculated for a given threshold on the 
coherent network amplitude rj. To display the depen- 
dence on component masses, the visible volume is binned 



8 


(25Mo x25M 0 bins) in the component mass plane. Here, 
we also assume that the detection efficiency is averaged 
over the sky position, binary inclination, and polarization 
angles. 

Instead of a direct calculation of Eq. 4.1, the integral 
can be estimated as a sum over the inverse density of the 
detected injections. Namely, each injection is assigned 
a density number pi and the analysis of each injected 
event is a statistical trial of whether or not that density 
is detected for a given threshold on 77 . The integral 4.1 
becomes then a sum over detected injections 

Kia (mi, m 2 ,r?) = ^ 4nrf (^r 1 ^,)) < 

(4.2) 

where r* is the distance to the i th injection, and dN^/dr 
is the radial density of simulated events. They are in- 
jected into a spherical volume with a fiducial radius of 
2 Gpc with the density distribution linearly increasing in 
distance and optimized to reduce the statistical errors. 
To express the search sensitivity, below we also . use the 
effective range which is calculated as the radius of 
the visible volume. 


C, False Alarm Rate Density and Event 
Significance 


The methods employed in previous searches for calcu- 
lation of event significance compare foreground events to 
the expected background. Given a foreground event with 
the coherent amplitude rj , its significance is determined 
by the false alarm rate 


FAR(ij) = 


N(y) 

2"bkg 


(4.3) 


where 7b kg is the accumulated live time for the corre- 
sponding background sample and N ( 77 ) is the number of 
background events with strength greater than 77 . How- 
ever, the IMBH search combines searches from two dif- 
ferent detector networks. Therefore, the coherent net- 
work amplitudes calculated for different networks are not 
directly comparable and the networks may have signifi- 
cantly different sensitivities and background rates. To 
combine the results of multiple searches into a single mea- 
surement, the IMBH search employs a statistical proce- 
dure based on the false alarm density (FAD) rate [54] 
defined as 


FAD ( 77 ) = 


1 

2"bkg 



(4.4) 


Given an event with the coherent amplitude 77 , its FAD 
rate is calculated from the mass averaged visible volume 
Vvis as a function of the coherent network amplitude. The 
sum is performed over the background events from the 
corresponding network with rji > rf. The FAD estimates 


the number of background events expected in a given net- 
work’s visible volume. Whereas the FAR statistic takes 
into account only the background rates, the FAD statis- 
tic also includes the sensitivity of a network to a popula- 
tion of expected GW sources [55, 56]. It weights search 
networks by their overall sensitivity to the source popu- 
lation and their background rates, and therefore allows a 
direct comparison between disparate networks. The FAD 
statistic is then used to rank candidates in the combined 
search. More significant candidates have smaller FAD 
rates. To determine the event significance, its FAD rate 
is compared to the time- volume product of the combined 
search (the overall search productivity v): 

l/(FAD) = ]TT obs [fc]VWfc](FAD), (4-5) 

k 

where the sum is over the networks and T 0 bs is the ob- 
servation time of each network (listed in Table I). The 
product p = FAD • z/(FAD) is the mean number of events 
expected from the background Poisson process. The false 
alarm probability (FAP) is calculated as 

N-l n 

FAP(A0 -l-££r exp(-M) . (4.6) 

n=0 

where N is the number of foreground events below a given 
FAD value. The FAP value indicates the probability that 
the candidate events are originating from a non-GW pro- 
cess. 


D. Statistical and Systematic Errors 

There are several uncertainties associated with the es- 
timation of the visible volume such as statistical errors 
due to a limited number of simulated events, calibration 
errors of the detector data streams, and systematic errors 
due to uncertainties of the simulated IMBH waveforms 
arising from differences between the waveforms and na- 
ture. 

By using binomial statistics and Eq. (4.2) the statis- 
tical uncertainty on the visible volume can be estimated 
as 

(4 - 7) 

where the sum is taken over the detected injections. The 
approximations used in the calculation could only in- 
crease the uncertainty, and therefore the estimate (4.7) 
is conservative. The statistical error in any given com- 
ponent mass bin is usually less than 5%. 

The calibration procedure of the GW strain data and 
the associated uncertainties for the S5/VSR1 run are de- 
scribed elsewhere [57, 58]. The amplitude calibration er- 
ror [ 20 , 21 ] directly translates into the error on the effec- 
tive range of the search (<11%). Respectively, the error 
on the visible volume is approximately 33%. 
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Further checks were performed using an updated 
EOBNR (EOBNRv2) [59] family which includes more 
PN corrections. Comparisons between the EOBNRv2 
family, the IMRPhenom family, and waveforms drawn 
directly from numerical relativity simulations agree to 
within 15% in the SNR induced in the detectors. Propa- 
gating this to the volume by noting that SNR is propor- 
tional to distance, we estimate a conservative systematic 
uncertainty on the search visible volume of 45% due to 
imperfect knowledge of the IMBH waveforms used in the 
simulations. 

The EOBNR waveform family (EOBNRvl) [12] used 
to calculate the visible volume in this search predicts 
more GW power radiated during the inspiral and merger 
phase than seen in numerical relativity simulations. As a 
result, a priori this model allows sources that are farther 
away to be seen. To account for the overestimation of the 
visible volume, the search for IMR signals in the S5 data 
between 25-100 M 0 [18] applied a distance correction. 
We follow the same procedure in the IMBH search. The 
newer EOBNRv2 family also includes higher order PN 
corrections to the dominant (2,2) mode and an improved 
calibration of the frequency evolution which manifest as a 
systematic shift in its noise weighted power as a function 
of frequency and a slight time dependence to the polar- 
ization. It is observed as an additional 10% bias in terms 
of 7] and a corresponding loss of detection efficiency due 
to the constant polarization constraint. For the case near 
the most sensitive region in the component mass space 
(centered around 88-1-88 M 0 ), the correction reduces the 
effective range of the search by about 50% . The effect on 
larger total mass systems is even more pronounced. 

It should be noted the EOBNR model used in the 
simulations neglects contributions to the waveform from 
other multipolar harmonics: Z,m modes different from 
the dominant 2,2 harmonic were not included in the ini- 
tial EOBNR model. These modes could provide addi- 
tionally detectable power. However, they interfere both 
constructively and destructively with the dominant har- 
monic and this may break the exact elliptical dependence 
(i.e. 90° phase shift) between the two polarization states. 
Overall, the effect on the detected power by these modes 
is expected be small in regards to other sources of uncer- 
tainty [60]. 

Propagating all the remaining uncertainties considered 
here into the volume gives an overall uncertainty of 60% 
in volume. The rate density estimates are then read- 
justed upward accordingly by the same amount. 


V. RESULTS AND DISCUSSION 

No event candidates were found to be significant to 
claim a detection. Therefore, we place upper limits on 
the rate density of IMBH coalescences as a function of 
the component masses. 


A. Event Candidates 

The FAD rate distributions for the background and 
foreground events are shown in figure 1 as a function of 
the coherent network amplitude. All events are ranked 
by their FAD rate, with the most significant events repre- 
sented by the low FAD values. Several foreground events 
with the lowest FAD rates are shown in Table IV. The 
first (rank 1) event with the lowest FAD rate is produced 
by the four detector network with rj = 3.16. This event 
has an associated FAP of 45%, which is not considered 
to be significant. No other event candidates have a low 
FAP sufficient for detection. 



FIG. 1: False alarm density rate vs 77 for the background 
events (H1H2L1V1 - solid line, H1H2L1 - dashed line) and 
the foreground events (H1H2L1V1 - black squares, H1H2L1 - 
open squares). 


B. Visible Volume and Rate Limits 

The visible volume is calculated from equation 4.2 for 
the events binned in the component mass plane. The 
thresholds on the coherent amplitude of each network 
are determined by the FAD rate of the loudest event (see 
Figure 1), denoted below in the text as FAD*. Figure 2 
shows the effective range as a function of the component 
masses for the networks analyzed in this search. The 
mass bins are limited to mass ratios less than 4:1, since 
no numerical relativity data is readily available for val- 
idation of the waveforms with larger mass ratios. For 
the more sensitive H1H2L1V1 network the best effec- 
tive range is achieved in the 88+88 M 0 bin at 241 Mpc. 
For the H1H2L1 network the corresponding range is 190 
Mpc. The ranges in Figure 2 take into account the SNR 
bias correction for the EOBNR waveform family as de- 
scribed in section IV D. Combining errors from the statis- 
tical procedure, calibration and the waveform systematic 




10 


TABLE IV: Highest ranked events by FAD. The first ranked event, produced by the four detector network, has a relatively 
small 7 compared to the other three events. However, the four detector network is much less noisy resulting in a low FAD 
value. 


Tank GPS time 

network 

V CC A„ et 

FAR (yr _ 

') FAD (Mpc - 3 Myr~ 

x ) FAP 

1 

871474393 H 1 H 2 L 1 V 1 3.16 0.90 0.17 

0.76 

0.09 

45% 

2 

857692870 

H1H2L1 

3.T4 0.74 0.13 

1.61 

0.26 

63% 

3 

846735754 

H 1 H 2 L 1 

3.69 0.76 0.13 

1.91 

0.30 

45% 

4 

820091022 

H1H2L1 

3.55 0.83 0.05 

2.90 

0.42 

51% 


errors, the total uncertainty on the effective ranges is es- 
timated to be 20%. 

In the absence of detection, we set upper limits on the 
rate of IMBH mergers at the 90% confidence level by 
using the loudest event statistic [61]: 

2 3 

fi90% = i'(FAL)*)' • (5,1) 

The ^(FAD*) is the time-volume productivity of the 
search calculated at the FAD rate of the first ranked 
event (the top event in Table IV). The rate density upper 
limits calculated in a binning of the component masses 
are presented in Figure 3. The upper limit for the com- 
bined search, averaged over all masses, is estimated to 
be 0.9 Mpc -3 Myr _1 . In the most sensitive bin, the rate 
limit is nearly an order of magnitude greater than for the 
overall search at 0,13 Mpc -3 Myr~ 1 . 

Since globular clusters are the most likely hosts of 
IMBHs, we convert our overall search upper limit into an 
astrophysical density of 3xl0 3 GC -1 Gyr -1 . This rate is 
still few orders of magnitude above the predictions for 
IMBH-IMBH rates in [41]. It should be noted, however, 
that the predicted astrophysical rates are very uncertain 
due to the lack of knowledge of the distribution and for- 
mation of IMBH sources. 


VI. SUMMARY 

In this paper we have presented a search for gravita- 
tional waves from IMBH coalescences in the mass range 
of 100-450 M 0 and mass ratios up to 4:1. The search was 
performed in the S5/VSR1 data collected with two dif- 
ferent network configurations: H1H2L1 and H1H2L1V1. 
For identification of potential GW candidates we used the 
Coherent WaveBurst algorithm with a polarization con- 
straint for the first time. To establish the significance of 
candidates from either search network, we combined their 
detection statistics into a single measurement by using 
the false alarm rate density statistic. No plausible GW 
candidates were identified. From this search, we place 
upper limits on the rate density of the IMBH binaries as 
a function of the component masses. In the most sensi- 
tive mass bin (centered at 88+88 M 0 ) the rate limit is 
0.13 Mpc _3 Myr _1 . When averaged over the mass plane 


the rate limit is 0.9 Mpc _3 Myr” 1 at the 90% confidence 
level. 

The sensitivity of the search was estimated by Monte 
Carlo simulations of detection efficiency using waveforms 
from the EOBNR family with component masses uni- 
formly distributed on the mass plane. The most dom- 
inant source of error in this analysis is the systematic 
uncertainty (45%) due to accuracy of simulated IMBH 
waveforms used for the estimation of the search visible 
volume. There are a few features of black hole coales- 
cence which were neglected in the simulation studies, for 
instance, the effect of spinning component black holes. 
However, un-modeled searches like CWB are sensitive 
to the energy emitted in gravitational waves regardless 
of details of the waveform evolution. While the effects 
of spinning component masses in the binary have not 
been examined in this study in detail, it is expected that 
their inclusion could only increase the effective range 
of the search [62], This is because the increase in the 
GW energy output in favorable (co-aligned) spin config- 
urations is greater than its decrease from less favorable 
(anti-aligned) configurations, giving an overall increase in 
the emitted energy. Moreover, the additional energy of 
aligned spin configurations could extend the mass range 
for which we can search beyond 450 Mq . 

This search has been limited to a relatively small area 
in the component mass plane of potential IMBH sources. 
However, future experiments with advanced detectors 
will have a significant increase in sensitivity, and more 
importantly, advanced detectors will also widen the sen- 
sitive frequency band. At low frequencies, LIGO design 
sensitivity at 10 Hz calls for an increase of a few orders 
of magnitude; hence a greater fraction of IMBH binary 
signals should become observable. These improvements 
should allow for better chances of detection of IMBH 
sources. 


Acknowled gment s 

The authors gratefully acknowledge the support of 
the United States National Science Foundation for the 
construction and operation of the LIGO Laboratory, 
the Science and Technology Facilities Council of the 
United Kingdom, the Max-Planck-Society and the State 
of Nieder sachsen/ Germany for support of the construe- 



11 




250 


200 


5 150 


100 


50 100 150 2UO 25U 3UU 35U 

Mass 1 (MJ 


200 


5 150 

CM 

(A 

<a 

* 100 


50 100 150 200 250 300 350 

Mass 1 (MJ 


- 200 

• j 180 
- 160 

140 

* - 120 


- H 100 

- 80 
■ 60 



FIG. 2: The effective range « in Mpc (values and white-gray scale) per component mass bin: H1H2L1V1 (left plot), H1H2L1 
(right plot). The overall uncertainty on the ranges quoted here is 20%. IMBH systems with mass ratios of greater than 4:1 are 
excluded since their waveforms are not verified with the NR calculations. 



200 


2 

CM 

(A 

tfi 

to 


50 


100 


200 250 300 350 

Mass 1 (M e ) 


1000 

100 

10 

1 

0.1 


FIG. 3: The rate density upper limits in Mpc“ 3 Myr _1 (values 
and white-gray scale) calculated over the component mass 
plane. These rate limits include all sources of error. 


tion and operation of the GEO 600 detector, and the Ital- 
ian Istituto Nazionale di Fisica Nucleare and the French 
Centre National de la Recherche Scientifique for the con- 
struction and operation of the Virgo detector. The au- 
thors also gratefully acknowledge the support of the re- 
search by these agencies and by the Australian Research 
Council; the Council of Scientific and Industrial Research 
of India, the Istituto Nazionale di Fisica Nucleare of 
Italy, the Spanish Ministerio de Education y Ciencia, 
the Conselleria d’Economia Hisenda i Innovacio of the 
Govern de les Illes Balears, the Foundation for Funda- 
mental Research on Matter supported by the Nether- 
lands Organisation for Scientific Research, the Polish 
Ministry of Science and Higher Education, the FOCUS 


Programme of Foundation for Polish Science, the Royal 
Society, the Scottish Funding Council, the Scottish Uni- 
versities Physics Alliance, the National Aeronautics and 
Space Administration, the Carnegie Trust, the Lever- 
hulme Trust, the David and Lucile Packard Foundation, 
the Research Corporation, and the Alfred P. Sloan Foun- 
dation. This document has been assigned LIGO Labora- 
tory document number PI 100068. 


Appendix A: Likelihood Analysis of Elliptically 
Polarized Waves 


The CWB algorithm performs the likelihood analy- 
sis [46] of the detector data streams, which are trans- 
formed into the time-frequency domain with the Meyer 
wavelet [63]. The sampled network data is x[z] = 
[xi[i], x/r[i]]j where k is an index iterating over the 
K detectors in the network. The vector x[z] is a func- 
tion of the time and frequency indicated with a single 
time-frequency index z, which is often omitted later in 
the text. The likelihood ratio is defined as 


A( x ,n) 


p(x|h(n)) 

p(x|0) 


(Al) 


where H is a parameter set describing two GW polariza- 
tions h = (h+, h x ) : p(x|0) is the joint probability that the 
data is only instrumental noise, and p(xjh) is the joint 
probability that a GW signal h is present in the data. 
The explicit form of the likelihood ratio is determined 
by the noise model p(x|0) and by the signal model h(D). 
In the analysis we assume that the noise of detectors is 
Gaussian with the standard deviations tr*. [i] . To account 
for the time and frequency variability of the noise, the 
cr fc [i] are estimated for every time-frequency sample. For 
elliptically polarized waves originating at a sky location 
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(0,<£), instead of reconstructing two unknown signal po- 
larizations A+[t] and (un-modeled case), only one 
waveform h[i] and two other signal parameters need to 
be reconstructed: the ellipticity parameter a (related to 
the inclination angle of the binary axis) and the polariza- 
tion angle Therefore the signal model is introduced 
into the anal^is by the following parameterization of the 
detector response 

€hW = *)M*1 + afx (fl, ®)*[i] , (A2) 

€h W = f+(n» *)*M - afx (n, > (A3) 

where A and £ h are the 90° phase shifted counterparts 
of h and £ h . The components of the response vector 
£ h [i], the noise scaled antenna pattern vectors f+( X )[i] 
and the noise scaled network data vector w[x] represent 
the individual detectors: 


€hM 

— (Chi [*]i — 

i£hJt[*]) , 

(A4) 

f +(x) 

f A+(x) 
” V *i[»1 1 

fK+M \ 
)' 

(A5) 

w[t] 

_ 

— 1 r J 

• » r f * 

(A6) 



°Kb]J 


The solutions for A, a, and are obtained by variation 
of the combined likelihood functional £(w|£ h )+£(w|£ h ), 
where w is the quadrature of w. The likelihood func- 
tional £(w|£ h ) is defined as twice the logarithm of A 

£( w l£h) = 2 3 4 5 6 7 8 9 ( w l£h) - (ChICh) » ( A7 ) 


the equations A2 and A3. Respectively the reconstructed 
null stream is n[t] = w[i] — £[i]. 

The three data streams w, £ and n are used to cal- 
culate the normalized (by the noise variance) energies of 
the network data stream E t otl the signal energy Eg w and 
the noise energy E n respectively. The Eq w = (£|£) = 
Ylij Liji where the components of the likelihood matrix 
Lij are calculated from the output of the corresponding 
detector pair (i, j). The sum of the off-diagonal terms 
(i jk j ) of the matrix Lij define the network coherent en- 
ergy E c = The phase shifted coherent energy 

E c is similarly defined for the phase shifted likelihood. 
There are two coherent statistics, that are used for CWB 
selection cuts: the network correlation coefficient 


cc = 


E c + E c 


(A9) 


\E C + E c | + E n + E n 
and the network energy disbalance A ne t — max(A,A), 
where 


A 



(A10) 


and E Ci E n: and A are calculated for the 90° phase shifted 
data streams. The CWB algorithm also defines the re- 
duced coherent energy 


e - V L V 


(All) 


where the inner products are given by a sum over a time- 
frequency area / containing the signal 

(a|b)=£(a[*]b[i]). (A8) 

*€/ 

The reconstructed network responses f [t] and f [i] are ob- 
tained by substituting the solutions for h, a and # into 


which is used for calculation of the main CWB detection 
statistic — the coherent network amplitude 

/W- < ai2) 

It is used to rank CWB events and establish their signif- 
icance against a sample of background events. 
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